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ABSTRACT 



Context. Recent studies of low- and intermediate-mass stars show that the evolution of the chemical elements in these stars is very 
different from that proposed by standard stellar models. Rotation-induced mixing modifies the internal chemical structure of main 
sequence stars, although its signatures are revealed only later in the evolution when the first dredge-up occurs. Thermohaline mixing 
is likely the dominating process that governs the photospheric composition of low-mass red giant branch stars and has been shown to 
drastically reduce the net 3 He production in these stars. The predictions of these new stellar models need to be tested against galaxy 
evolution. In particular, the resulting evolution of the light elements D, 3 He and 4 He should be compared with their primordial val- 
ues inferred from the Wilkinson Microwave Anisotropy Probe data and with the abundances derived from observations of different 
Galactic regions. 

Aims. We study the effects of thermohaline mixing and rotation-induced mixing on the evolution of the light elements in the Milky 
Way. 

Methods. We compute Galactic evolutionary models including new yields from stellar models computed with thermohaline instability 
and rotation-induced mixing. We discuss the effects of these important physical processes acting in stars on the evolution of the light 
elements D, 3 He, and 4 He in the Galaxy. 

Results. Galactic chemical evolution models computed with stellar yields including thermohaline mixing and rotation fit better ob- 
servations of 3 He and 4 He in the Galaxy than models computed with standard stellar yields. 

Conclusions. The inclusion of thermohaline mixing in stellar models provides a solution to the long-standing " 3 He problem" on a 
Galactic scale. Stellar models including rotation-induced mixing and thermohaline instability reproduce also the observations of D 
and 4 He. 

Key words. Galaxy : evolution - Galaxy : abundances - Galaxy : formation. 



1. Introduction 

Understanding the evolution of the light elements deuterium (D), 
helium-3 ( 3 He), and helium-4 ( 4 He) hinges on the comprehen- 
sion of several astrophysical processes and links together dif- 
ferent branches of physics and cosmology. D, 3 He, and 4 He 
are all synthesized in astrophysic ally relevant quantities during 
Big B ang Nucleosynthesis (BBN; |Peeblesi ri966; Wa goner et alj 
I 967|). In the absence of any other realistic pro duction channel 
( Epstein et al.lll976tlProdanovic & Fieldsll2003l) . the abundance 
of D in galaxies smoothly decreases in time as gas cycles through 
stars. 3 He and 4 He, instead, have a more complex history, since 
they are both produced and destroyed in stars. 

The sensitivity of the predicted BBN abundance of D, 
(D/H)p, to the baryon-to-photon ratio, t/b, coupled with its 
straightforward galactic evolution, has long made th is ele- 
ment be the 'baryometer' of choice dReeves et al.lll973l) . In the 
nineties, the controversial assessment of the primordial abun- 
dance of D from observations of high-r edshift clouds (e.g. 
ISongaila et al. I ll994t iBurles & Tvtler|[l998l) . probing an almost 



unevolved medium, led researchers to resort to detailed Galactic 
chemical evolution (GCE) modeling in order to constrain the 
pre-Galactic D abundance and, hence, the value of It was 
shown that any reasonable GCE model r esults in moderate local 
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local present-day abundance of D is well known, this provided 
a stringent bound to (D/H)p and, hence, a test for BBN theories. 
Chemical evolution models able to reproduce the majority of the 
observational constraints for the Milky Way f irmly ruled out the 
highest values of (D/H)p by ISongaila et al.l (1 19941). In par ticu- 
lar, the values of (D/H )p and 7/b suggested byljosi (2000J) and 
IChiappini et al.l (120021) turned out to be in very good agreement 
with the ones determined from the analysis of the first Wilkinson 
Microwave Anisotr opy Probe (WMAP) data dBennett et alj2003fc 
LSoereel et al.ll2003h . 

Indeed, the first release of results from WMAP made us 
enter a new precision era for cosmology - t]b is now known 
with exquisite accuracy. Converging measurements of D abun- 
dances in_jejriote_gas_clouds lead to (D/H)p = (2.8 + 0.2) x 
10~ 5 dPettini et alj 120081) . a value consistent, within the er- 
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Table 1. Abundances of Deuterium, 3 He, and 4 He at different epochs. 



Chemicals 


Units 


SBBN+VKMAP 
(t=0 Gyr) 


Protosolar cloud 
(t=9.2 Gyrs ) 


Local intestellar medium 
(Present time) 


D 


10 5 (D/H) 


2.49 ±0.17" 

2 60 +0A9b 

z " u " -0.17 


2.1 +0.5 C 


2.31± 0.24" 
0.98 ±0.19'' 
2.0 ±0.1^ 


J He 


10 3 ( j He/H) 


1.00 ±0.07° 
1.04 ± 0.04* 


1.5 ± 0.2 C 


2.4 ±0.7« 


4 He 


Y 


0.2486 ± 0.0002° 

f) 7A7Q +0 00046 


0.2703" 





(a)ICyburt et al] d2008h: (b)ICoc et all (l2004l>: <c) iGeiss & Gloeckle j lll998h : (d)lLinskv et al] 
j20061) ; (e)lHebrard et al.l d2005h ; (fl lProdanovic et aljfeOlCft : (g) lGloeckler & Geissldl996l) ; 
(h) lAsplund et all l l2009h 



rors, with the primordial D abundance predicted by the stan- 
dard mo del of cosmology with parameters fixed by WMAP 
data (e.g. lCvburt etafll2008l: ICoc et al.ll2004 see Table Q). The 
remarkable homogeneity of D abundances at high redshifts, 
however, clashes with the unexpectedly large scatter in D/H 



(Vidal-Madiaretal. 


1998;Uenkinsetal. 1999b 


Sonneborn et al] 


2000: Hebrard et al. 


20021 Oliveira & Hebrard 


2006k. The ob- 



served dispersion can be reconciled with the predictions on D 
evolution from standard GCE models by taking into account 
two short-term, small-scale phenomena: D deple tion onto dust 
grains and localized infall of primordia l gas dRomano et al.l 
120061 ISteigman et aljl2007t lR oman oll20Toh . 

As far as 3 He is concerned, liber] d 19671) and 
iTruran & Cameronl (1 1 97 lb first showed that large amounts 
of this element are produced by low-mass stars (M 1-3 
M ) in the ashes of hydrogen burning by the p- p cycle on 
the m ain sequence. Problems arised soon, when iRood et all 
dl976l) incorporated these yields in models of GCE and came 
out with predicted 3 He abundances orders of magnitude higher 
than the observed ones in the Galaxy. Indeed, the nearly 
constancy of the 3 He ab undance with both time and position 
within the Galaxy (e.g. iBania et ail 120021) rather implies a 
negligible production of this element in stars, at variance with 
predictions from standard stellar models. It was then advocated 
that some non-standard mechanisms are acting in a major 
fraction of low-mass stars, which prevents the fresh 3 He from 
surviving and be ing ejected in the interstellar medium (ISM; 
IRood et al.l[l984h . For instance, extra mixing was also invoked 
in order to explain other abu ndance anomalies l i ke the carbon 
isotopic ratio in red giants (ICharbpnnell Il995l Hoganl 19951: 
Charbonnel & Do Nascimentol Il998t see also lEggleton et al.l 
2006L 120081) ." The first 3 He stellar yields for low-mass stars 
computed with "ad hoc" extra-mixing (i.e., not related to 
a physical mecha nism, see Section 2.2.2) became available 
with the work of iBoothrovd & Sackmannf (119991) . iGalli et al.l 
d!997l) : iPalla et al] d2000l) : IChiappini et al] (l2002h implemented 
those yields, accounting for the dependency of the helium-3 
production/destruction on the stellar mass, metallicity and initial 
D abundances. It was shown that chemical evolution models 
which account for about 90% of low-mass stars undergoing 
extra-mixing led to a good agreement with the Proto Solar 
Cloud (PSC) observations as well as with the observed gradient 
along the disk. 

The physical process responsi ble for the extra-mixing has 
possibly been recently identified. ICharbonnel & Zahnl (l2007bl) 
have shown indeed that thermohaline instability, when mod- 
elled with a simple prescription based on linear stability anal- 



ysis, leads to a drastic reduction of 3 He production in low- 
mass, low-metallicity red giant stars. However a couple of plan- 
etary nebulae, namely N GC 3242 and J32 0, have been found 
to behave classically (see lBania et al EoiOh : slightly more mas- 
sive than the Sun, they are currently returnin g fresh 3 He to the 
ISM, in ag reement with standa rd pre dictions (IRood et aDl 19921 : 
IGalli et all 119971: iBalser et al.l 1 19991 120061) . To reconcile the 
3 He/H measurements in Galactic HII regions with the high val- 
ues of 3 He in NGC 3242 and J320. ICharbonnel & Zahnl d2007al) 
proposed that thermohaline mixing is inhibited by a fossil mag- 
netic field in red giant branch (RGB) stars that are descendants 
of Ap stars. The percentage of such stars is about 2-10% of all 
A-type objects. This number agrees with the fact that about 4% 
of low-mass evolved stars (including NGC 3242) exhibit stan- 
dard sufa£e_abundanc£s^ carbon isotopic ra- 
tios d Charbonnel & Do Nascimentolll998l) . Grids of 3 He yields 
from stellar models taking into account thermoh aline instabil- 
ity, as well as rotational mixing, are presented in ILagarde et al.l 
d201 ll) . These are suitable for use in GCE studies such as the one 
presented here. 

The scarce determinations of reliable 4 He abu ndances in the 
Galaxy dBalser et al.ll2010t iPeimbert et~ai1l2010l) . together with 
the steadiness of GCE model predictions on the evolution of 4 He 
when adopting differ ent stellar yields (e.g. IChiappini et al.l2002l: 
iRomano et al.ll2.Q10l) . have given this element little attention in 
GCE studies. 

In this paper, we deal with the evolution of D, 3 He, and 
4 He in the solar vicinity, as well as their distributions across 
the Galactic disc. The production of 3 He is strictly related to 
the destruction of D and the production of 4 He. Therefore, these 
elements are considered all together. The main novelty of the 
present work is the new nucleosynthesis prescriptions for the 
synthesis of 3 He and 4 He in low-mass stars (masses below 6 
Mo). Here we present, for the first time, chemical evolution 
models for 3 He computed with stellar yields from non-standard 
stellar mode ls that include both rotation-induced and thermoha- 
line mixing dCharbonnel & Zahnll2007aHCharbonnel & Lagardel 
12010b ILagarde et al.l 1201 lb . As we will show in the next sec- 
tions, these yields bring chemical evolution model predictions 
into agreement with the Galactic 3 He data without the necessity 
of assuming "ad hoc" fractions of extra-mixing among low-mass 
stars. 

The layout of the paper is as follows. In Section [2] we dis- 
cuss the production/destruction channels of D, 3 He, and 4 He in 
stars in light of new generation stellar models that take into ac- 
count rotation-induced mixing and thermohaline instability. A 
comparison with stellar yields predictions from the literature is 
presented. In Section [3] we describe the adopted GCE model. In 
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Section|4]we discuss our results on the evolution of the light el- 
ements in the Galaxy. In particular, we try to put constraints on 
the uncertain parameters of stellar evolution through a compari- 
son of the model predictions with the relevant data. We draw our 
conclusions in Section[5] 



2. Stellar nucleosynthesis of light elements from 
new generation models for low- and 
intermediate-mass stars 



In a series of paper s (ICharbonnel & Lagardel 120101 
lLagarde et all 1201 ll 120121, hereinafter Papers I II and III, 
respectively; see also Charbonnel & Zahr3l2007blah . we discuss 
the impact of rotation-induced mixing and thermohaline insta- 
bility on the structure, evolution, nucleosynthesis and yields, as 
well as on the asteroseismic and chemical properties of low- 
and intermediate-mass stars at various metallicities. A detailed 
description of the input physics of the models is given in Paper 
III. Rotation-induced mix ing is treate d usin g the complete 
forma lism developed by IZahnl d 19921) and iMaeder & Zahnl 
dl998l) (see for more details, Papers I, II, III). In these stellar 
models, we consider that thermohaline instability develops 
as long thin fingers with the aspect ratio consistent with 
predictions b y lUlrichl (1 1972b and confirmed by the laboratory 
experiments dKrishnamurtil 120031) . We note that this value is 
higher t han obtained by cu r rent 2D and 3D numerical simu- 
lations (|Denissenkovl 120101; Denissenkoy & Merrvfie'ldl 1201 It 



iRosenblum et al.ll201 lUTraxler et al.ll201 l}. Before a final word 
on this discrepancy comes from future numerical simulations 
in realistic stellar conditions, we adopt an aspect ration (i.e., 
maximum length relative to their diameter) of 5, which nicely 
accounts for the observed chemical properties of red giant stars. 
These new generation stellar models account very nicely for 
main-sequence and RGB abundance patterns observed in field 
and open cluster stars over the m ass and metallicity range cov - 
ered, as shown in Paper I and in ICharbonnel & Zahnl (l2007bb . 
In this section we briefly summarize their characteristics as far 
as the nucleosynthesis of light elements is concerned. 

2.1. Deuterium 

Deuterium burns by proto n-captures at low t emperature dur- 
ing the pre-main sequence dReeves et al. I fl97l . Therefore, it is 
totally destroyed in stellar interiors and the corresponding net 
yields are negative whatever the mass and metallicity of the star, 
as in the case of standard stellar models. 



2.2. Helium-3 

2.2.1 . Impact of rotation-induced mixing and thermohaline 
instability on model predictions 

Paper II describes in detail the behaviour of 3 He both in the 
standard case and in our models including thermohaline insta- 
bility and rotation-induced mixing. Figure Q] shows the theo- 
retical evolution of the 3 He surface abundance along the evo- 
lutionary tracks in the HR diagram for the models of 1.5, 3.0, 
and 6.0 M at solar metallicity that take into account these two 
processes. The colour coding in mass fraction is given on the 
right, with the initial values of 3 He and D+ 3 He indicated by the 
arrows. Evolution is shown from the pre-main sequence along 
the Hayashi track up to the end of the thermally-pulsing asymp- 
totic giant branch (TP-AGB) phase. We see the changes in 3 He 



surface abundance due to D-burning on the pre-main sequence 
except in the 6.0 M model whose convective envelope with- 
draws very quickly at the beginning of that phase allowing to 
the preservation of pristine D in a very thin external layer. 3 He 
is then produced in low- and intermediate-mass stars during the 
main sequence through the pp-chains and subsequently dredged- 
up when the stars move towards the RGB. On the 1.5 M track 
one sees the effect of rotation-induced mixing that already brings 
fresh 3 He towards the stellar surface while the star is on the main 
sequence (on the contrary, in the standard case the surface abun- 
dance of 3 He remains constant on the main sequence and starts 
changing only at the very base of the RGB due to the first dredge- 
up). However, as discussed in Paper II rotation-induced mixing 
is found to lower the total 3 He production compared to the stan- 
dard case over the whole mass and metallicity range scrutinized, 
and to decrease the upper mass limit at which stars destroy this 
element (see Fig. 5 in Paper II). Additionally, for low-mass stars 
(M < 2-2.2 M ) thermohaline mixing occurring during the RGB 
phase beyond the bump and on the TP-AGB leads to the de- 
struction of part of the freshly produced 3 He while accounting 
for the observed surface abundance anomalies of other chemi- 
cals (i.e ., lithium and nitrogen, as w ell as carbon isotopic ratio; 
see also ICharbonnel & Zahn 2007b); the associated decrease of 
surface 3 He is clearly seen for the 1.5 M model in FigQ] (right 
panel). Therefore, although low-mass stars remain net 3 He pro- 
ducers, their contribution to the Galactic evolution of this ele- 
ment is much lower than in the standard case. It was also shown 
in Paper II that thermohaline mixing leads to 3 He depletion dur- 
ing the TP-AGB phase for stars with masses up to ~ 4 M (see 
the 3 M track in the figure). In more massive intermediate-mass 
stars, 3 He is further destroyed through hot-bottom burning on 
the TP-AGB. Finally, this figure shows that the lower the mass 
of the star, the higher the surface abundance of 3 He at the end 
of the TP-AGB. The global impact of rotation-induced mixing 
and thermohaline instability on the net 3 He yields from low- and 
intermediate-mass stars of various metallicities is summarized in 
Fig. 9 and in Tables 1 to 4 of Paper II (see also Fig.|2]discussed 
below). 

2.2.2. Comparison with other stellar models 

In Fig. [2] we show the mass fraction of 3 He at the end of the 
second dredge-up as a function of the initial stellar mass at two 
metallicities (Z=0.0001 and 0.014; left and right panels, respec- 
tively) from our models (black lines; see Paper II). A compari- 
son is made with model predictions from the literature (coloured 
lines and filled triangles). 

Our standard predictions (full lines) are in very 



(1996). 

IBoothrovd & Sackmannl(ll999l) . and lWeiss et all (119961) . For the 



good agreement with those fro m iDearborn et al 



reasons given in Section 12.2.11 they are higher than in the case 
including rotation-induced mixing and thermohaline instability 
(dashed black lines). Let us note that thermohaline mixing on 
the TP-AGB leads to further decrease of the 3 He mass fraction 
at the stellar surface (reduction of 64, 83, and 19% in the 
[1.25 M , Z ], [0.85 M , Z=0.0001], and [2.0 M , Z and 
Z=0.0001] models, respectively; see Paper II). 

In Fig. [2] we also plot the predi ctions at the end of the 
second dredge-up for the models of IBoothrovd & Sackmannl 
(I1999L hereafter BS99) that include parametric mixing below 
the base of the convective envelope of RGB stars (the so-called 
"conveyor-belt" circulation and the associated "cool bottom pro- 
cessing", CBP; dashed magenta lines). In BS99 post-processing 
computations, the mixing is not related to any physical mech- 
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log Teff (K) log Teff (K) 

Fig. 1. Theoretical evolutionary tracks in the HR diagram for the 1.5 M Q , 3.0 M , and 6.0 M models at solar metallicity following 
the standard prescription (left panel) ; and including rotation-induced mixing and thermohaline instability (right panel), from the 
pre-main sequence up to the end of the TP-AGB phase. Colours depict the mass fraction of 3 He at the stellar surface as indicated on 
the right, with the arrows showing the initial 3 He and 3 He+ 2 H content assumed at birth. 



anism. Rather, the depth of the mixed zone is a free parame- 
ter that corresponds to the difference AlogT between the tem- 
perature at the base of the hydrogen-burning shell and that at 
the base of the assumed mixed zone. BS99 assumed that the 
value of AlogT remains constant along the RGB, and this free 
parameter was calibrated in order to reproduce the observed car- 
bon isotopic ratio of clump stars in M67. We note that the fixed 
AlogT value used by BS99 leads to shallower mixing just after 
the RGB bump than what we get when considering thermoha- 
line instability as the physical mixing mechanism. As a conse- 
quence BS99 prescription leads to a slow and gradual theoret- 
ical decrease of the 12 C/ 13 C ratio up to the RGB tip; this be- 
haviour disagrees with the sudden drop of this quantity that is 
observed after the RGB bump and that is well reproduced by 
our models that include the rmohaline mixing (see Paper I and 
ICharbonnel & Zahdl2007bl However, due to the higher com- 
pactness of the hydrogen-burning shell when the stars reach 
higher luminosities on the RGB, the fixed AlogT value adopted 
by BS99 leads to deeper mixing close to the RGB tip. The sec- 
ond free parameter is the stream mass flow rate, M p . Contrary to 
thermohaline mixing, for which the stream mass flow rated de- 
creases along the RGB from ~10~ 2 5 to ~10~ 5 M Q yr -1 , BS99 
assumed that the value of M p stays constant along the RGB 
(M p = lO _4 M yr _1 ). This difference in the mixing efficiency 
at the end of the RGB explains the difference between BS99's 



1 To translate the thermohaline diffusion coefficient of our models 
into a stream mass flow rate as used by BS99, we compute : 



with r and p the radius and the density at the base of the region where 
the thermohaline instability develops, /„„,- the size of this mixing zone, 
and D the mean diffusion coefficient along the mixing zone. 



models and ours on the final 3 He surface mass fraction. Another 
drawback of such models was also that a physical reason was 
missing to explain why only 7% of the low mass stars would 
produce 3 He following the standard predictions. 

2.3. Helium-4 

2.3.1. Standard predictions 

Changes in the surface 4 He abundance result from the dredge- 
up episodes undergone by the stars along their evolution, that 
eventually lead to the partial engulfment of the ashes of central- 
and shell-H-burning by the convective envelopes. The efficiency 
of the successive dredge-up events depends both on stellar mass 
and metallicity (see e.g. Paper II), as depicted in Fig.[3]where we 
present our standard predictions (upper panels) for the surface 
abundance of 4 He, adjusted to its initial value, in mass fraction, 
after the first and the second dredge-up (solid and dashed lines, 
respectively) as a function of the initial stellar mass, for our two 
extreme metallicities (Z and Z=0.0001). For solar metallicity, 
the surface abundance of 4 He increases during the first dredge- 
up over the whole mass range investigated, while for the lowest 
metallicity (Z=10~ 4 ) only stars with M < 2.5 M are affected. On 
the other hand, the second dredge-up leads to an increase of the 
4 He surface abundance only for the upper part of the considered 
mass range (i.e., for stars with masses higher than ~ 2.5 M ). 
Finally, the third dredge-up, which occurs in stars with mass > 
6.0 M forZ and 2.5 M forZ=10~ 4 during the TP-AGB phase, 
allows to increase the 4 He abundance at the surface of these stars 
(not shown here). This increase is ~ 1% for [6.0 M , Z ], and 
[4.0 Mq, Z=0.0001] modelsQ 



2 Note that we did not include additional mechanisms to force the 
third dredge-up in our models, as required in th e literature for re- 
producing the carbon-star luminosity function (e.g. iFrogel et al.lll99d : 
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Z = 0.0001 Z=0.014 




M(M G ) M(M S ) 

Fig. 2. Mass fraction of 3 He at the stellar surface at the end of the second dredge-up as a function of initial stellar mass in standard 
models (solid lines and triangles) and in models including various prescriptions for mixing in the radiative regions (dashed lines), 
at two metallicitie s (Z=0.0001 and Z=0.014; left and right panels, respectively). The predictions are from iDearborn et al.l (l!996[ 
DS T96. blue line) iBoofhrovd & Sackmannl(ll999l BS99, at the RGB tip, magenta lines). ! Weiss et alJ(ll996l WWD96, red triangles), 
and lLagarde et al.l (1201 ll LCI 1. black lines). The initial abundances of 3 He adopted at the zero age main sequence by the different 
authors are given on the top right corner of each panel. 



Z e , stand. 4jj e Z e 4 He 




4.2 4 3.8 3.6 3.4 4.2 4 3.8 3.6 3.4 

log Teff (K) log Teff (K) 



Fig. 5. Same as Fig.Q]for 4 He from the zero age main sequence to the end of the TP-AGB. 



ICosta & Frogellll996UGroenewegenlll998l;lMarigo et al.lll999h . Hence, 
we do not attempt any comparison of the model predictions with avail- 
able data for C and N isotopes. However this has no impact on the galac- 
tic chemical evolution of the light elements we focus on. 



2.3.2. Impact of rotation-induced mixing and thermohaline 
instability on model predictions 

Thermohaline instability has no significant impact on the sur- 
face abundance of 4 He since it develops only in the upper wing 
of the hydrogen-burning shell (see e.g. Paper I for more de- 
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Fig. 3. Relative enrichment of the surface abundance of 4 He with 
respect to its initial value (in mass fraction) for Paper I models at 
the end of the first and second dredge-up (solid and dashed black 
lines, respectively) as a function of the initial stellar mass and 
for two metallicities (Z and Z=0.0001 ; left and right panels, re- 
spectively). Top panels: our standard predictions are compared 
with those of BS99 (blue circles; open and full at the end of the 
first and second dredge-up, respectively). Bottom panels: our 
models including both rotation-induced mixing and thermoha- 
line instability are compared with the "cool-bottom processing" 
predictions by BS99 at the RGB tip (blue squares). 

tailsB However, rotation-induced mixing smoothes the inter- 
nal abundance gradients compared to the standard case, leading 
freshly produced 4 He to diffuse outwards. This can be seen in 
Fig. |4] where we show the abundance profile of 4 He at the end 
of the main sequence for two stellar masses and metallicities in 
the standard and the rotating cases (black full and red dashed 
lines, respectively; the vertical bars indicate the depth reached 
by the convective envelope at its maximum extent during the first 
dredge-up). This results in stronger 4 He abundance variations at 
the end of the first dredge-up when rotation-induced mixing is 
accounted for, as seen in Fig. [3] (bottom panels). For the same 
reasons, the effect of the second dredge-up is also strengthened 
in the rotating models. Overall, the impact of rotation increases 
with decreasing metallicity (see Paper I, and references therein). 

Consequently, the yields of 4 He are higher in the models in- 
cluding rotation-induced mixing than in the standard case (see 
tables [3] to [6). And at the same time, the fraction of H is de- 
creased. 

Figure [5] summarizes all the effects described above on the 
surface abundance of 4 He along the evolution of standard and 
rotating stars with different initial masses at solar metallicity. 
The tracks are shown from the zero age main sequence up to the 
end of the TP-AGB, with the colour coding given on the right 
scale. 

3 Thermohaline mixing leads to an increase of the surface abundance 
of 4 He during the RGB by ~0.1% in the 1.25 M Q models independently 
of the stellar metallicity. 



Fig. 4. Abundance profile of 4 He at the end of the main sequence 
for the 1 .5 M and 2.0 M models at two metallicities (Z=0.0001 
and Z=0.014, left and right panel, respectively) in the standard 
case (solid black line), and when including the effects of rotation 
(red dashed line). Models are from Paper II. The maximum depth 
reached by the convective envelope during the first dredge-up is 
shown with the vertical bars and arrows. 

2.3.3. Comparison with other stellar models 

In Fig. [3] we compare ou r predictions with those of 
iBoothrovd & Sackmannl d 1999b . We note a very good agreement 
as far as the standard models are concerned (upper panels). On 
the other hand, neither the thermohaline instability nor BS99 
"cool bottom processing" do affect the 4 He surface abundances. 
The differences that appear in the lower panels result essentially 
from the effects of rotation-induced mixing. 

3. Chemical evolution model 

3.1. Basic assumptions 

The adopted model for the chemical evolution of the Galaxy 
assumes that the Milky Way forms out of two main accretion 
episo des almost completely disentangled dChiappini et al.llT997l 
During the first one, the primordial gas collapses very 
quickly and forms the spheroidal components, halo and bulge. 
During the second one, the thin-disc forms, mainly by accre- 
tion of matter of primordial chemical composition. The disc is 
built-up in the framework of the inside-out scenario of Galaxy 
formation, whi ch ensures the formation of abundance grad ients 
along the disc dLarsonll 19761 iMatteucci & Francoislll989l) . The 
Galactic disc is approximated by several independent rings, 2 
kpc wide, without exchange of matter between them. 

3.2. Nucleosynthesis prescriptions 

The nucleosynthesis prescriptions for metals (elements heavier 
than 4 He) are from Paper II for low- and intermediate-mass stars 
(M < 6 M ). As for massive stars and Type la supernovae, we 
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Tabl e 2. Description of d ifferent models computed with yields 
from lLagarde et all (1201 ll) for stars below 6M . 



Models 


Stellar mass 




M< 2.5 M Q 


M>2.5 M Q 


A 


100 % standard 


100 % standard 


B 


96 % th. + rot llJ 


100 % th.+ rot 




4 % standard 




C 


100% th.+ rot 


100 % th.+ rot 



th. + rot means that we use yields including 
thermohaline instability and rotation induced 
mixing. 



adopt the same nucleosynthesis prescriptions as in lRomano et al.l 
(120101) . their model 6, namely: 

- Yields for cor e- collap se supernovae are taken from 
iKobavashi et alj ( 120061) . except for carbon, nitrogen 
and oxygen, for which t he adopted y ields ar e from 
IMevnet & Maederl d2002l): IHirschi etail d2005l) : iHirschil 
d2007l) : lEkstrom et al.l d2008l) . 

- Yields for Type la supernovae are from llwamoto et al.l 
(fl999h . 

As far as the nucleosynthesis prescriptions for D, 3 He 
and 4 He are concern ed, we use standard prescriptions from 
iDearborn et al.l (1 1996b in the 6 to 100 Mr mas s range for D 
and 3 H e and yields by IMevnet & Maedej d2002l ): IHirschi et alj 
d2005l) : IHirschil (12007k : lEkstrom et al l (l2008l) for 4 He. For low- 
and intermediate-mass stars (M < 6 M ) we use the yields from 
Papers II and III0 and we compute three different models with 
the following assumptions (summarized in Tabled : 

- In Model A we consider the standard yields for all stars with 
masses below 6.0 M Q . 

- In Model B we consider the yields including thermoha- 
line instability and rotation-induced mixing for 100% of the 
intermediate-mass stars (M > 2.5 M ) and 96% of the low- 
mass stars (M < 2.5 M ). The remaining 4% of the low-mass 
stars are assumed to release the yields as predicted by the 
standard theory. This allows us to account for the fact that 
some rare planetary nebulae (J320, NGC 3242) exhibit 3 He 
as pr edicted by the standard models (see e.g. iBalser et"a71 
120071 and references therein). 

- In Model C we assume the yields including thermohaline 
instability and rotation-induced mixing for all the low- and 
intermediate-mass stars. Obviously the outcome of models 
B and C are expected to be very similar. 

Linear interpolations in mass and metal licity are used to fill 
the gaps in the computed grids of yields (see lRomano et al.l2010l 
for a discussion of the potential spurious effects introduced by 
this procedure). 



4. Evolution of the light elements in the Milky Way 

In the following sections we discuss the evolution of the light 
elements in the Galaxy within the framework described in 
Section [3] when taking into account the yields from our non- 
standard stellar models, and compare the predictions with dif- 
ferent observations. The relevant data are presented in TableQ] 



4. 1 . Evolution of deuterium 

In this paper we assume (D/H)p = 2.6 x 10~ 5 for consistency 
with the initial value adopted by the stellar models. However, 
we note that our GCE model gives a better fit to the data when 
a value of (D/H)p - 2.8 x 10 as observed in quasa r spectra 
dPettini et al.ll20 08) is adopted instead (see Fig. 3 of iRomanol 
2010). 

Figure [6] (left panel) shows the data for deuterium from the 
Big Bang (WMAP value, filled black diamond) to the present 
day (local interstellar medium, LISM, filled triangles). The 
determinations of prim ordial and PS C (filled green circle, 
Geiss & Gloecklerl fl998l according with iGeiss & Reeveslfl972l 



19811) deuterium abundances underline a small depletion from 
the Big Bang (t=0 Gyr) to the solar-system formation (t=9.2 
Gyr). Analyses of Far Ultraviolet Spectroscopic Explorer 
(FUSE) observations have allowed measurements of D/H in 
the LISM for many lines of sight. These observations have 
revealed a large variation of the local D abundances, which 
complicates t he interpretation i n the contex t of standard 
GCE mode ls dLinskv et alJl2006t lLinskvll2010l) . IHebrard et all 
d2005l) and lLinskv et all d2006T) ~proposed that either the lowest 
(D/H=0.98 x 10" 5 ) or the highest (D/H=2.31 x 10~ 5 ) observed 
value is indicative of the true LISM value reflecting the process 
of D astration through succe ssive stellar g e neratio ns during 



r generatioi 
~aD d200l 



4 3 He stellar yields are presented in Table 1 to 4 of Paper II. 4 He 
stellar yields are shown in Tables [3]to[6] 



the whole Galaxy's evolution. IHebrard et al.l (12005D suggest a 
value of the true local deuterium abundance lower than the one 
measu red in the local bubble. On the other hand, iLinskv et al.l 
J2006h give a lower bound to the true local deuterium abundance 
very close to the primordial abundance, pointing to a deuterium 
astration factor smaller than pred i cted b y standard GCE models. 
More recently, IProdanovic et al.l (1201 Ol) have applyed a statisti- 
cal Bayesian method to determine the true local D abundance. 
They propose a value very close to the D abundance at the time 
of the formation of the Sun (see Table [TJ. GCE models that 
fulfil all the major observational constraints available for the 
solar neighbourhood and for the Milky Way disc can explai n the 
LISM D abundance suggested by IProdanovic et al.l d2010h as a 
result of D astration during Galactic evolution; lower and higher 
values can only be explained as due to small-scale, transient 
phenomena, such as D depletion on to dust grain s and localized 
infall of gas of primord ial chemical com positio n dRomano et al.l 
I2006h ISteigman etail 120071: IRomanol 1201 Ol and references 
therein). 

Our predictions for the evolution of deuterium with time 
in the solar neighborhood and for the present deuterium abun- 
dance profile along the Milky Way disc are shown in Fig. [6] 
(left and right panels, respectively). The very modest shift be- 
tween the predictions of Models A and B is due to the slight 
difference in the H yields when stellar rotation is accounted for 
compared to the standard case. Since pristine deuterium burns 
entirely in stars independently of their mass and metallicity, the 
predicted trends depend only on the total astration assumed in 
the GCE model. The present predictions fit well the PSC deu- 
terium da ta and favour the local abundance of deuterium sug- 
gested by IProdanovic et all (l2010l) . 



4.2. Evolution of He 

Figure [7] (left panel) shows the evolution of 3 He/H in the so- 
lar neighborhood as predicted by Models A, B, and C (see 
Sect. 13.21 1. As already discussed in the literature (see Sect. [TJ, 



7 



N. Lagarde et al.: Evolution of light elements in the Galaxy : D, 3 He and 4 He 



T 



~i — i — i — r 



T 



~i — i — i — r 



Solar neighbourhood 



J_ 



j i i i_ 



5 10 
Time (Gyr) 



15 




10 

R G (kpc) 



Fig. 6. Evolution of D/H in the solar neighborhood (left panel) and distribution of D/H along the Galactic disc at the present time 
(right panel). The predictions of Models A and B are shown (red solid and dash ed blue lines, respectively). At t=0 Gyr, we plot the 
WMAP value (filled black dia mond). The PSC data bvlGeiss & Gloecklerl d!998|) are shown (filled gree n circle). The local interstellar 
medium (LISM) data are from lLinskv et al.l d2006l) . lHe brard et al.1 (I2005T) and lProdanovic et aD ( 20 1 0|) (filled orange, light blue, and 
magenta triangles, respectively). The data for the outer disc (Rc=10 kpc) are from lRogers et al. 1 (120051 filled bordeaux circle). 



standard stellar models strongly overestimate the production of 
3 He by low-mass stars. 

According to standard predictions, 3 He is produced on the 
main sequence during the core hydrogen burning and it is not 
destroyed during later phases. Consequently, when adopting 
standard predictions 3 He is overproduced in the course of 
Galactic evolution (Model A, red solid line). In low-mass 
stellar models including thermohaline mixing, 3 He is destroyed 
from the bump luminosity on the RGB during shell hydrogen 
burning and during the early AGB. When low-mass stars 
are assumed to experience this extra-mixing, GCE models 
do not overproduce 3 He (Fig. UJ Model C, cyan dotted line). 
Model B (dashed blue line) shows the effect of inhibiting the 
thermohaline mixing in 4% of low-mass stars (see Section EOl ). 
that thus follow the standard prescriptions. As can be seen 
in Fig. |7J negligible differences are found between Model 
B and Model C predictions. The PSC and local 3 He a bun- 
dances dGeiss & Gloecklerl 1 19981 iGloeckler & Geissl 1 19961 ) are 
reproduced by the models (at 3- and 1-sigma level, respectively). 

The present-day 3 He abundance distribution along the 
Galactic disk is shown in Fig. |7j (right panel). Model A, adopt- 
ing standard 3 He yields for low-mass stars, predicts too much 
3 He at all Galactocentric distances. In the inner regions, where 
the star formation was stronger, the contribution of low- and 
intermediate-mass stars to the 3 He enrichment of the ISM is 
more important. As a consequence, Model A predicts a large 
negative gradient of 3 He/H. Although the contribution of low- 
mass stars in the inner regions is significantly reduced in Model 
B, it still predicts a negative gradient. In Fig. |7j we com- 
pare the predictions of Models A and B wi th observations of 



3 He in HII regions fro m iBania et alj (2002) and with the lo- 
cal 3 He abundance from IGloeckler & Geissl ( 19961) . Contrary to 



what is predicted by classical theory, t he observat i ons in HII 
regions show a gradient close to zero. Ba nia et all d2002b de- 
rived meaningful error estimates only for one source [S209 : 
3 He/H=(l.l+0.2) x 10~ 5 ) and pointed our the riskiness of bas- 
ing one's conclusion on only one object. In order to better con- 
strain both stellar and galactic evolution studies, it would be 
very useful if future work could provide a sound estimate of 
the errors associated to 3 He determinations across the Galactic 
disk. Nevertheless, we can conclude that thermohaline and ro- 
tating mixing are required to fit currently available measure- 
ments of 3 He in Galactic HII regions. In addition, the 3 He gra- 
dient is dlog( 3 He/H)/dR c 0.04 and -0.028 dex/kpc for 

mo dels A and B, r espect ively, consistent with that predicted 
bv IChiappini etal] d2002h . -0.04<dlog( 3 He/H)/dR G < -0.03 
dex/kpc. 

iTosil (120001) showed that G CE models including the CBP 
dBoothrovd & Sackmannl 19*991) also fit the relevant Galactic 3 He 
data, pro vided the primordial abundance of D is suffi ciently low 
(see also IChiappini et a"fll2002t iRomano et alJl2003h . However, 
CBP destructs 3 He too quickly and is not a physical mechanism 
(see Section l2~2l i. Therefore, it does not provide a reliable phys- 
ical explanation of the 3 He observations in HII regions. 

4.3. Evolution of 4 He 

Figure [8] shows the temporal (left panel) and the spatial (right 
panel) evolution of the helium isotopic ratio ( 3 He/ 4 He) in the 
solar neighborhood and along the Galactic disk at the present 
time, respectively. Only the predictions of Models A and B 
are s hown in the right panel. The PSC and LISM data ar e 
from iGeiss & Gloecklerl d 19981) and IGloeckler & Geissl (Il996l) . 
respectively. The standard theory predicts a significant evolution 
of the helium isotopic ratio in the Galaxy, and a strong gradient 
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Fig. 7. Left panel : Evolution of 3 He/H with time in the solar neighborhood. Data for the PSC (green fil l ed cir cle) and local 
interstellar medium (LISM, orange filled triangle) are from lGeiss & Gloecklerl ( 1 19981) and iGloeckler & Geissl (1 19961) . respectively. 
1-0" and 3-cr error bars are shown with thi ck and thin l i nes, re spectively. Right panel : Radial distribution of 3 He/H at the present 
time. The dots are HII regions data from lBania et al. (120021) (error bars are shown only for S209; see text for discussion). The 
triangle at Rg=8 kpc represents LISM data from IGloeckler & Geissl (Il996t . The predictions from Models A, B and C are shown in 
both panels by the red full, blue dashed and cyan dotted lines respectively. 



across the Galactic disc. When the thermohaline mixing - which 
destroys 3 He in the stellar interior - and the rotation-induced 
mixing - which increases 4 He at the stellar surface (see Section 
3.3) - are taken into account, lower helium isotopic ratios are ob- 
tained (cfr. the dashed versus solid lines in Fig. [8}. Consequently, 
the gradient of the helium isotopic ratio across the Galactic disc 
decreases and the observations in HII regions can be reproduced. 

Figure [9] presents the predictions for Y versus [Fe/H] in the 
solar neighborhood. The model predictions are comp ared with 
the in itial abundance of 4 He in the Sun derived by ISerenellil 
(120101) from full standard solar models computed for differ- 
ent sets of solar abundanc es. When the solar abundances from 
iGrevesse & Sauval (1 19981) and a metal to hydrogen ratio in the 
solar photosphere of (Z/X) p i,-0.0229 are adopted, a value of 
Y,-=0.2721 is obtained, which results in good agreement with 
the predictions of chemical evolution models including rotation- 
induced mixing and thermohaline mixing (green Sun symbol 
and dashed and dotted lines in FigO resp ectively). On the other 
hand, choosing the solar abundances from lAsplund et aU (120091) 
and (Z/X) p/l =0.0178 result in ¥,=0.2653, which is better fit- 
ted by models including standard stellar nucleosynthesis (orange 
Sun symbol and solid line in FigO. Since the actual solar chem- 
ical composition is matter of debate, we can not discriminate 
between our models basing on the results for Y. 



5. Conclusions 

In this article, we have described the results for the Galactic 
evolution of the primordial elements D, 3 He and 4 He in the 
light of new stellar models including the effects of thermohaline 
instability and rotation-induced mixing on the stellar structure, 
evolution and nucleosynthesis. The new stellar models are 



presen ted and discussed at len gth in ICharbonnel & Lagardel 
d2010h and lLasarde etal.1 d2011l) . We model the time-behavior 
of D, 3 He and 4 He in the solar neighbourhood, as well as 
in the inner and outer disc. The predictions of our models 
including thermohaline instability and rotation-induced mixing 
are compared to the predictions of models adopting standard 
nucleosynthesis prescriptions and to the most recent relevant 
observations. 



We have seen that the CBP proposed before by 
iBoothrovd & Sackmannl (1 19991) is not a physical mechanism, 
and yields deduced do not take into acc ount phases beyond 
the RGB. As discussed in the literatu r e (ICharbonnel & Zahnl 



l2007bt ICharbonnel & Lagardel l2010t lLagarde et al.l 120111) 



thermohaline mixing can induces significant depletion of He in 
low- and intermediate-mass stars. Although these stars remain 
net producers of 3 He, their contribution to the Galactic evolution 
of this element is highly reduced compared to classical theory. 
Indeed, our GCE models including thermohaline mixing repro- 
duce the observations of 3 He in the PSC, LISM and HII regions, 
while 3 He is overproduced on a Galactic scale with standard 
models. Thermohaline mixing is the only physical mechanism 
known so far able to solve the so-called " 3 He problem" plaguing 
in the literature since many years. Importantly, its inhibition by 
a fossil magnetic field in red giant stars that are descendants of 
Ap stars does reconcile the measurements of 3 He/H in Galactic 
HII regions with high values of 3 He in a couple of planetary 
nebulae. 

On the other hand, rotation has an impact on the stellar yields 
of H and 4 He. Whether GCE models including rotation fit 
better the Galactic 4 He data than standard models can not be 
said, because of the uncertainties on the actual solar chemical 
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Fig. 8. Temporal {left panel) and spatial {right panel) variation of helium isotopic ratio in the solar neighborhood a nd along the 
galactic disk at the present time, respec tively. Models are the same as in Fig|7] Dat a are from Geiss & Gloecklerld 19981) for the PSC 
(green filled circle). iBania et al.1 (|2002|) for Galactic HII regions (black dots), and lGloeckler & Geissl (1 1996b for the LISM (orange 
filled trangle). 1-cr and 3-<x errors are also shown as thick and thin lines, respectively. 
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Fig. 9. Y versus [Fe/H] in the solar neighborhood predicted 
by models including yields computed with standard prescrip- 
tions (Model A; solid line) and with thermohaline mixing and 
rotation-induced mixing (Models B and C; dashed and dotted 
lines, respectively). Mode l predictions are compared with the 
protosolar values given bv lSerenelhl (12010b . 



composition. However, they are co nsistent with the r elative ly 
high value of (D/H) L ism proposed bv lProdanovic et all (1201 Oh . 



We conclude that GCE models including both thermohaline 
mixing and rotation-induced mixing reproduce satisfactorily 
well all the available data on D, 3 He and, possibly, 4 He 
abundances in the Milky Way, within the errors. However, 
the additional data on HII regions and better refining of their 
interpolation would be crucial to assess how good our fit 
actually is. 



Acknowledgements. We dedicate this paper to Robert T. Rood : the scientist to 
whom the stellar evolution and 3 He communities owe so much, and the friend 
whom we will miss for ever. We thank T. Bania, D.Balser, and H.Reeves for 
fruitful discussions. We are thankful to our referee, Dr Achim Weiss, for his 
interesting and constructive remarks on our article. C.C., D.R., and M.T. grate- 
fully acknowledge the enlightening and fascinating conversations with Johannes 
Geiss and his hospitality at the International Space Science Institute (ISSI) in 
Bern (CH). We acknowledge financial support from the Swiss National Science 
Foundation (FNS) and the french Programme National de Physique Stellaire 
(PNPS) of CNRS/INSU. 



References 

Asplund, M„ Grevesse, N, Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Balser, D., Rood, T. R., & Bania, T. M. 2010, in IAU Symposium, Vol. 268, IAU 

Symposium, ed. C. Charbonnel, M. Tosi, F. Primas, & C. Chiappini, 101-106 
Balser, D. S., Bania, T. M., Rood, R. T„ & Wilson, T. L. 1999, ApJ, 510, 759 
Balser, D. S., Goss, W. M., Bania, T. M„ & Rood, R. T. 2006, ApJ, 640, 360 
Balser, D. S., Rood, R. T., & Bania, T. M. 2007, Science, 317, 1171 
Bania, T. M., Rood, R. T., & Balser, D. S. 2002, Nature, 415, 54 
Bania, T. M., Rood, R. T., & Balser, D. S. 2010, in IAU Symposium, Vol. 268, 

IAU Symposium, ed. C. Charbonnel, M. Tosi, F. Primas, & C. Chiappini, 

81-90 

Bennett, C. L., Halpem, M., Hinshaw, G., et al. 2003, ApJS, 148, 1 
Boothroyd, A. I. & Sackmann, I.-J. 1999, ApJ, 510, 232 
Buries, S. & Tytler, D. 1998, ApJ, 507, 732 
Charbonnel, C. 1995, ApJ, 453, L41 + 

Charbonnel, C. & Do Nascimento, Jr., J. D. 1998, A&A, 336, 915 
Charbonnel, C. & Lagarde, N. 2010, A&A, 522, A 10+ 
Charbonnel, C. & Zahn, J. 2007a, A&A, 476, L29 



10 



N. Lagarde et al.: Evolution of light elements in the Galaxy : D, 3 He and 4 He 



Table 3. Model results for metallicity of Z=0.0001 
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Each row contains entries for different assumptions: standard 
(without fhermohaline or rotation-induced mixing); thermoha- 
line mixing only; thermohaline and rotation-induced mixing. 

1 The initial rotation on the ZAMS 

2 Yields of 4 He 

Table 4. Same as Table[3]forZ=0.002 
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1.90. 10~ U2 




111. T^IUL. 


U.H- J 




J.JJ. 1U 


2.5 


stand. 






2.94. 10"" 2 




th. +rot. 


0.45 


141 


5.62. 10- 02 


3.0 


stand. 






2.03. 10^ U2 




th. +rot. 


0.45 


147 


6.75. 10-° 2 


4.0 


stand. 






7.83. 10"" 2 




th. +rot. 


0.45 


147 


1.93. 10-°' 


6.0 


stand. 






3.49.10-'" 




th. +rot. 


0.45 


167 


4.14.10- 01 



Table 6. Same as Table[3]for Z=0.014 



M 




VzAMsA'cr/r 


VzAMS 


Yield 4 He 


(M G ) 






(km.sec -1 ) 


(M ) 


1.0 


stand. 






1.29. 10-" 2 


1.25 


th. +rot. 


0.45 


110 


2.27.10-" 2 


1.5 


stand. 






2.79 .10-" 2 




th. +rot. 


0.45 


110 


2.45 .lO^ 02 


2.0 


stand. 






1.65. 10-" 2 




th. +rot. 


0.45 


110 


3.36. 10-° 2 


2.5 


stand. 






2.35. 10-" 2 




th. +rot. 


0.45 


130 


5.36. 10- 02 


3.0 


stand. 






4.08. 10-" 2 




th. +rot. 


0.45 


136 


8.09. 10-° 2 


4.0 


stand. 






4.80. 10-" 2 




th. +rot. 


0.45 


144 


8.98. 10- 02 


6.0 


stand. 






3.07.10-'" 




th. +rot. 


0.45 


156 


3.84.10-'" 
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